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INTRODUCTION 


The  purpose  of  this  report  is  to  give  a brief  account  and  summary 
of  the  research  accomplishments  under  the  present  AFOSR  Grant  to 
Texas  A & M Research  Foundation.  The  unsteady  airload  prediction 
technique  and  an  aeroelastic  analysis  of  a cascade  of  staggered 
blades  in  subsonic  compressible  flow  developed  under  this  research 
grant  are  discussed.  Several  research  papers  and  reports  have  been 
published  from  the  research  conducted  and  a list  of  these  papers 
and  the  summaries  are  included  in  the  Appendix. 

Many  researchers  have  devoted  considerable  effort  towards  predicting 

the  flow  characteristics  through  the  multiple  stages  of  blades  that 

exist  in  axial  flow  compressors.  The  research  efforts  have  been 

concentrated  on  obtaining  the  aerodynamic  loads  utilized  in  the 

design  of  more  efficient  blades.  In  addition  to  the  need  for 

information  about  blade  aerodynamic  loading,  flutter  characteristics 

also  need  to  be  determined.  With  the  assumption  of  a two-dimensional 

flow,  the  flow  through  a single  row  of  blades  is  mathematically 

equivalent  to  the  flow  through  a staggered  cascade  of  infinitely 

1 2 

many  airfoils.  Early  researchers,  Whitehead  , Kemp  and  Sears  , 

3 

and  Schorr  and  Reddy  , assumed  incompressible  flow  and  developed 
aerodynamic  theories  and  computational  schemes  for  predicting  the 
unsteady  airloads  on  the  blades  for  oscillatory  freestream  flow 

2 

and/or  oscillating  blades.  An  example  is  the  work  of  Kemp  and  Sears 
who  studied  the  problem  of  the  unsteady  lift  generated  on  a reference 
airfoil  of  a cascade.  In  their  approach,  they  used  an  oscillatory 
freestream  flow.  Their  study  considered  the  steady  interaction 
between  blades  but  neglected  the  unsteady  interaction  and  hence, 
the  effect  of  cascade  spacing.  Their  approach  was  to  express  the 
unsteady  lift  as  a function  of  the  design  parameters,  such  as  the 
ratio  of  the  airfoil  chord  and  the  disturbance  wavelength,  thus 
enabling  a designer  to  optimize  the  performance  of  a turbomachine 
design  instead  of  analyzing  a particular  blade  arrangement.  In 

3 

another  study,  Schorr  and  Reddy  treated  the  flow  through  a stag- 
gered cascade  of  airfoils  in  which  the  effect  of  unsteady  upstream 
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disturbances  were  included  as  in  the  case  of  unsteady  or  distorted 

inlet  flow  conditions  in  an  axial  flow  compressor.  Their  problem 

was  also  formulated  under  the  assumption  of  incompressible  potential 

flow  and  numerical  results  were  obtained  for  oscillatory  flow  using  an 

approximate  solution  developed  from  the  integral  equations  involved. 

in  addition,  their  solution  yielded  unsteady  lift  coefficients  for  the 

airfoils  as  a function  of  the  frequency  of  the  oscillations  and  for 

different  values  of  stagger  and  solidity  of  the  cascade. 

4 

In  an  independent  study,  Jones  and  Moore  studied  the  incompressible 

flow  about  a cascade  of  oscillating  airfoils  at  zero  mean  incidence. 

In  order  to  obtain  a solution,  they  utilized  a unique  numerical 

lifting  surface  technique  which  differs  from  all  other  methods  in 

that  it  makes  use  of  the  velocity  potential  instead  of  acceleration 

5 

potential  doublet  distributions.  Rao  and  Jones  later  applied  this 

technique  to  the  oscillatory  flow  about  an  airfoil  of  a staggered 

cascade.  Airload  results  obtained  for  several  values  of  frequency, 

interblade  spacing  and  stagger  angle,  showed  excellent  agreement 

3 

with  the  results  of  Schorr  and  Reddy  . 

Adopting  a technique  similar  to  that  of  Schorr  and  Reddy,  Fleeter^ 
considered  the  effects  of  compressibility  on  both  the  fluctuating 
lift  and  moment  coefficients  for  cascaded  airfoils  having  an  upstream 
non-uniformity.  He  obtained  a solution  for  the  time-dependent,  two- 
dimensional,  partial  differential  equations  which  describe  the  per- 
turbation velocity  potential  through  an  application  of  Fourier 
transform  theory.  The  resulting  integral  solution  equation  was 
evaluated  numerically  by  a matrix  inversion  technique.  The  fluctuating 
lift  and  moment  coefficient  variations  were  computed  and  represented 
as  a function  of  Mach  number,  cascade  solidity,  cascade  stagger 

angle,  interblade  phase  lag,  and  reduced  frequency.  Jones  and 

7 .... 

Moore  extended  the  velocity  potential  formulation  to  oscillating 

two-dimensional  airfoils  in  compressible  flow.  In  their  numerical 

method,  they  replaced  the  slowly  converging  Hankel  function  series 

by  a rapidly  converging  exponential  series.  They  studied  the  effects 

of  varying  airfoil  spacing,  frequency,  Mach  number,  and  phase 


difference  between  adjacent  blades.  Variations  in  the  aerodynamic 

damping  can  become  zero  but  never  negative  at  certain  discrete 

frequencies.  This  is  a desirable  characteristic  with  respect  to 

flutter  due  to  bending.  The  results  also  indicated  that  the 

pitching  moment  aerodynamic  damping  can  become  zero  relative  to 

the  blade  quarter-chord  axis,  while  also  being  zero  at  the  critical 

frequencies,  and  could  be  negative  at  the  higher  Mach  numbers  over  a 

wide  range  of  frequencies  of  interest  in  flutter  analysis.  This 

is  an  undesirable  characteristic  from  the  standpoint  that  it 

increases  the  area  of  instability  for  torsional  flutter.  Rao  and 
0 

Jones  utilized  the  theory  developed  in  Ref.  7 to  determine  the 
airload  and  moment  coefficients  on  a reference  airfoil  of  a stag- 
gered cascade  of  airfoils  in  subsonic  flow.  Circumferential  dis- 
torsion  due  to  inflow  conditions  was  expressed  as  an  interblade  phase 
lag.  Results  were  obtained  for  several  values  of  frequency,  Mach 
number,  interblade  spacing,  stagger,  and  phase  lag  angles  for  both 
cases  of  oscillatory  flows  and  oscillating  blades.  The  oscillatory 
flow  results  compared  well  with  those  of  Fleeter. 

9 

In  his  study  utilizing  compressible  flow,  Whitehead  presented 
calculations  for  the  torsional  flutter  of  a cascade  of  unstalled 
blades  at  zero  mean  deflection  and  subsonic  Mach  numbers.  Whitehead 
found  that  the  effect  of  increasing  Mach  number  was  favorable  and 
tended  to  suppress  the  flutter  that  was  predicted  by  incompressible 
theory . 

In  this  report  the  general  aerodynamic  theory  and  the  numerical 
lifting  surface  theory  is  presented.  An  investigation  is  conducted 
for  a single  degree  of  freedom  system  in  torsion.  The  effect  of 
the  flow  and  geometric  parameters  is  evaluated  in  establishing 
flutter  boundaries  and  these  results  are  compared  with  those  of 

9 

Whitehead  who  used  a completely  different  computational  procedure 
for  calculating  the  aerodynamic  derivatives.  Additionally,  a 
general  flutter  program  is  applied  for  a two-degree  of  freedom 
(bending-torsion)  staggered  cascade  in  subsonic  flow.  By  utilizing 
an  iterative  procedure  which  permits  frequency  variation,  the  flutter 
frequency  and  the  flutter  speed  of  the  reference  airfoil  are  obtained 
as  a function  of  the  cascade  parameters. 
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AERODYNAMIC  THEORY 


General 


The  governing  equation  for  the  unsteady,  compressible,  two- 
dimensional  flow  of  an  isentropic,  inviscid,  irrotational  fluid 
is  given  in  terms  of  its  velocity  potential  by 


V2<(> 


l_r0  + |E^+5.v(§i)] 

cl  ot 


(i) 


where  <t>  is  the  perturbation  velocity  and 


<P+u)  1 + w jc  . 


(2) 


U is  the  freestream  velocity.  The  respective  perturbation  velocity 
components  along  the  x a z axes  are,  u (=|^-)  and  w • 

Assuming  £hat  u and  w are  small  compared  to  U,  Eqs . (1)  and  (2) 

are  combined  to  yield, 


(i-M2 ) + = (3) 

ax2  az2  a2  at2  a 8xat 

where  M is  the  freestream  Mach  number  and  a is  the  speed  of  sound. 
The  following  non-dimensional  coordinates  and  definitions  are  used: 

x = i'z  = ei'T  = ] r (4) 

<|>(x,z,t)  = U£<Kx,z)ei(eX+u/r^  , (5) 

where 

“ = Jp  t e = , and  # = 4>(X,Z)  . (6) 

U B2 
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When  Eq. (3)  is  combined  with  Eq.(5),  it  reduces  to  a two-dimensional 
Helmholtz  equation  for  the  pertubation  velocity  in  the  transformed 
coordinate  system, 

V2<J>  + k 2<i>  = 0 (7) 


where  k = . 

$2 

For  a flow  problem,  the  boundary  conditions  are  usually  prescibed . 

For  an  isolated  airfoil  in  compressible,  unsteady  flow,  the  solution 

to  Eq. (7)  can  be  derived  by  the  application  of  Green's  theorem. 

A relation  for  the  velocity  potential  at  a given  point,  $ , is 

P 

given  in  terms  of  velocity  potential  distribution  over  the  lift- 
ing surface  and  its  wake.  Treating  the  lifting  surface  as  a thin 
airfoil,  the  discontinuity  in  the  velocity  potential  between  the 
upper  and  lower  surfaces  is  expressed  as  a doublet  distribution, 

K(=  <J>  - ).  One  such  solution  for  an  isolated  airfoil 

upper  lower 

is  given  in  Ref. 7.  The  relation  between  the  downwash  at  any  point, 
p,  on  a thin  reference  blade  on  a staggered  cascade  in  subsonic 
flow  and  the  modified  doublet  distribution  K(X),  is  given  as 


2™p  = 


CO 

= -J  K(x] 


- S (X  -X,  Z , D,  H,  a ) dx 
, 2 u P P 


where, 


im (a+eD) 


X+mD)2  + (mH-zp)2 


]}  <9! 


j AW 

and  D = j,  H = -j-,  and  o is  the  interblade  phase  lag.  The  blades 
of  the  cascade  are  numbered  m,  with  m = 0 as  the  reference  airfoil. 
Since  SQ  in  Eq.  (9)  satisfies  the  wave  equation  in  the  form. 


+ ic2S_  = 0 
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and  in  the  limit  as  z -+0 , Eq.  (10)  may  be  rewritten  in  the  form, 

Jr 


2irWp  = / K (X) 


+ k S 1 dX 


where. 


;im(a+cD)  (Xp~x'*'rnD) 

(Xp-X+mD)  2+m2H2  ** 


(Xp-X+mD) 2 + (mH-D) 2 


The  series  involving  Hankel  functions  (SQ,S^)  in  Eqs.(9)  and 
(12)  have  very  poor  convergence  characteristics.  Therefore, 
these  are  replaced  by  an  exponential  series  as  shown  in  Refs.  10 
and  11.  However,  it  is  important  to  understand  that  this  trans- 
formation from  Hankel  function  series  to  exponential  series  is 
valid  only  for  an  infinite  cascade.  Hence,  it  cannot  be  applied 
to  a finite  cascade.  The  convergence  of  the  exponential  series 
is  so  good,  that  the  required  computational  time  is  less  for  a 
cascade  when  compared  to  a two-dimensional  isolated  airfoil  in 
subsonic  flow  where  it  is  required  to  use  Hankel  functions.  The 
transformed  relations  as  given  in  Ref. 7 are, 


so  - -5  E 


. oo  -27ra  (m)  | X -X  | / S 

£ £ t- 

2 m=-°°  (6-m)  2 _u 2 ** 


and. 


S1  " 


a(m)e-27Ta(m)|Xp-Xl/S 

~2—J  (6-m)  2-y2  ** 


, for  Xp  > or  < X , 


where , 


a (m)  = [ (6-m) *-y 2 ] ^ | + i(6-m)|  , for  Xp 


> or  < 
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and. 


6 


a+eD 
2u  ' y 


kS 

2it 


and  S 


(D^H2)*5  . 


(16) 


The  exponential  form  of  the  series  Sq  and  not  only  converge 
rapidly,  but  also  provide  directly  the  values  of  critical  frequencies 
as  shown  in  Ref. 7.  Eqs.(13)  and  (14)  will  diverge  whenever  one 

of  the  denominators  vanishes  and  no  solution  to  Eq. (11)  would  be 
possible.  The  critical  values  of  the  parameter  y for  which  the 
analysis  fails  is  given  by, 

V = 6,  1+6,  2+6,  etc. 

This  phenomenon  corresponds  to  a resonance  condition  at  critical 
frequencies  which  constitutes  an  infinite  set  of  values  of  a para- 
meter depending  on  flow  and  configuration  characteristics  at  which 

the  aerodynamic  function  becomes  infinite  everywhere.  Resonance 
conditions  are  functions  of  Mach  number,  frequency,  interblade 
spacing,  stagger  angle,  and  acoustic  velocity.  The;*  represent 
the  condition  at  which  self-induced  aerodynamic  forces  are  zero 
and  the  blades  act  effectively  as  if  they  were  in  a vacuum. 

Boundary  Conditions 

The  downwash  w(=  w,e1U)T)  can  be  expressed  in  terms  of  (X,Z)  by 
using  Eq. (5) , 

w = ||  = 6Uei  (tX+«.T)  u _ (17) 

Downwash  can  now  be  non-dimensionalized  to  give, 


W' 


8<J>  _ w,e-lex 
8z  BU 


(18) 


For  the  oscillating  blades  with  flapping  and  pitching  motions 
about  the  mid-chord  position,  the  downwash  boundary  condition 
is  defined  as, 


w'  ± = U [icoz ' + (l+iujxi ) a ' ) 


(19) 
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where  z'  and  a'  are  the  amplitudes  in  flapping  and  pitching  res 
pectively.  Since  periodic  motions  were  assumed,  z and  a are 
defined  as 


z 


z'eipt  = z ,eia)T 


a = a ' eipt  = a'ei(1)T 


The  doublet  distributions,  K 's,  are  complex  and  depend  on  z' 
and  ex'  for  given  values  of  w and  M.  A typical  will  be  of 
the  form 


Kn  " 


anz ' 


+ b a 1 
n 


where  a and  b are  complex  quantities  and  depend  on  frequency, 
n n 

Mach  number,  interblade  spacing,  stagger,  and  phase  lag. 

Numerical  Procedure 

The  numerical  procedure  developed  in  Ref.  7 is  utlized.  The  wake 
boundary  condition  is  given  by, 


K (X)  = K.  e_iv(X_1)  , 
te  ' 

where  K is  the  doublet  strength  at  the  trailing  edge  (X=l) . 
When  Eqs.  (11)  and  (23)  are  combined,  they  can  be  expressed  in 
the  form. 


where , 


dX  + KfceI 


I 


fe- iv(X-l) 

1 


K 


te 


r!!i 

[ax 


dX 


p 


ivSQt  - v2 (1-M2 )P 

Se-2Tra  (m)  (1-Xp)/S 
(2ira  (m) +ivs]  [ (6-m)  2-p2]Js 


r 


; 


i 


and  a(m)  = |(6  -m)  2 - vi2]1^2  ^ “ i(6-m)5.  The  symbols  SQt  and  Slfc 

represent  S and  respectively  at  the  trailing  edge. 

The  airfoil  is  divided  into  N equal  strips  and  K is  assumed  to  be 
Kn,  a constant  over  the  nth  strip,  where  n = 1,2, 3..., N.  Let  2B  be 
the  width  of  each  strip  and  denote  the  center  of  the  nth  strip. 
In  Eq.(24),  and  W are  replaced  by  X^  and  respectively,  where 
i refers  to  the  ith  strip.  Hence  Eq.(24)  will  then  be  given  by, 

N 

2*Wi  K"  1S1  <Xi’Xn-B>  ' S1(WB>  + 2*’BSo<Xi-Vl 

n=l 


+ Kte  I (Xi) , where  i = 1,2,..,N. 


(27) 


When  X.  = X , S is  evaluated  by  integrating  the  exponential  form 

of  S over  the  interval  Xn-B  £ x .£  Xn+B*  For  ot^er  intervals' 
X.^x°  and  S can  either  be  found  by  integration  or  be  taken  as  the 
mean^alue  of  SQ  over  the  interval.  Kfce  can  be  expressed  as  a 

function  of  Kv,  using  the  wake  condition.  This  relation  is  given 
N 


by, 


K 


K 


n 


te  (2ivB+e"lvB) 


(28) 


For  a given  geometry  and  flow  conditions,  the  numerical  terms  in 
Eq. (27)  can  be  evaluated  by  using  Eqs.(13),  (14),  (25),  and  (26). 

For  a known  set  of  W± ' s , Eq.  (27)  reduces  to  a system  of  N equations 

with  N unknowns,  1^,  K2,  K^,  when  it  is  combined  with  Eq. (28) . 

It  is  assumed  that  is  known  over  the  airfoil.  In  Eq.  (27)  , this 
represents  a set  of  linear  algebraic  equations  given  by , 


2ir{Wi}  = [A]  (K)  . 


(29) 


Therefore,  knowing  the  values  of  { wJ  and  [a]  allows  for  the 
determination  of  the  doublet  distribution,  (k). 
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Aerodynamic  Derivatives 

Euler's  equation  of  motion  in  unsteady  flow  is  given  as, 

|H  + u|u  = _1  l£  (3 

3t  u3x  p 3x  * 

In  terms  of  the  upper  and  lower  velocity  potential  and  pressure 

on  the  airfoil,  Eq.  (30)  can  be  shown  to  be, 

9k  . TT9k  1,  . , q 

3t  + U3?  = "p(p£“pu)  (3 

where  k(x)  = • Since  the  lift  per  unit  chord  per  unit  span 

is  £(x)  = (p  -p„)  , Eq.  (31)  becomes. 


£ ru' 


Now  from  Eq. (5) , 


e<!t  + “ft'  - *<*>  ■ 


k(x)  = = U£(4>u-4>£)el(eX+a)T) 


k (X)  = U£K(X)e1 (eX+wT) 


(33a) 


(33b) 


where  K(X)  = • Equations  (32)  and  (33)  are  combined  to 

yield, 


« , r 3K (X)  1 i (eX+wT) 

fc(X)  = pU2  JivK(X)  + -gj — ]e 

where  v=  u>/g2  Eq.  (34)  is  valid  on  the  airfoil  surface,  -1<  X <1. 
In  the  wake  region,  no  pressure  discontinuities  are  allowed  to 

exist  and  hence, 

1VK  (X)  + — — = 0 
must  be  satisfied  when  X = l. 

% 

When  values  of  K have  been  obtained,  the  local  lift,  £ (x)  * at  a 

^ _ i e x 

point  X is  given  by  Eq. (34) . If  K is  substituted  for  Ke  , then 


pU2  £ 


ft(X) 
pU2  £ 


t ‘ H + ll)et“T 


- I ifa>K  + 


(36a) 


(36b) 


The  lift  L(=  L'elpt)  and  the  nose-up  pitching  moment  about  the 
mid-chord  axis  M(=  M'elpt)  are  given  by, 


and, 


pU2£ 


- - Kte  i KdX  = C^z'  + Cta«- 


J.  -L 

= K + f KdX  - iio  f XKdX  = C z ' + C a' 

,2^2  te  J J mz  ma 


where  C . , C „ , C , and  C are  the  aerodynamic  derivatives . 

£z  Ha  mz  ma 

These  are  usually  complex  and  depend  on  the  geometry  and  flow  char- 
acteristics . 


EQUATIONS  OF  MOTION 

The  equations  of  motion  for  a two  degree-of-f reedom  (bending 
torsion)  system  of  an  airfoil  are 

M£z  + Sot  + Kz£z  = -pU2  i (C^zz  + C£aa) 


where. 


S«z  + 15  + Koc  = pu!*’<cm2z  + Cmaa) 


M =f dm  = blade  mass  per  unit  span, 

S = / r dm  = blade  static  moment  about  elastic  axis, 

*1  2 

I = I r dm  = blade  mass  moment  of  inertia  about  elastic 
J axis, 

A = semi-chord  length  of  the  blade, 
z = flapping  motion, 
a = pitching  motion, 

K • bending  stiffness  of  the  blade 
KQ  » torsional  stiffness  of  the  blade. 
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Assuming  simple 

harmonic  motion. 

2 = z’eipt  = z'e1" 

(41a) 

a = a'e^pt  = a,e^'a)T 

(41b) 

and  separating 

the  aerodynamic  derivatives  into  real  and 

imaginary 

parts , 

Clz  = Clzr  + lC£zI 

(42a) 

CA,a  = CSLaR  + lC£.aI 

(42b) 

Cmz  = CmzR  + lCmzI 

(42c) 

n — p x ^ C 

ma  maR  mal 

(42d) 

Eqs. (39)  and  (40)  reduce  to 


UKZ  - M&p2  + pU2JtC£z)z’  + (pU2  £C£a-Sp2 ) a ' = 0 
-(PU2£2Cmz  + SJtp2)z’  + (Ka  - pU2p2Cma-Ip2)a'  = 0 . 


(43) 


Equations  (43)  are  a set  of  two  linear,  homogenous  equations  in 
z'  and  a'.  For  this  system  to  have  a nontrivial  solution,  the 
coefficient  determinant  of  of  z'.and  must  be  equal  to  zero. 

UKz  - M£p2  + pU2i,C£z)  (pU2JtC£a  - Sp2) 

"(Cmzpu2£2  + S£p2)  (Ka  ' pu2*2cma  " Ip2) 

This  determinant  is  referred  to  as  the  flutter  determinant  and 
the  solution  gives  the  flutter  frequency.  However,  a direct 
solution  cannot  be  obtained  since  the  arodynamic  derivatives 
are  functions  of  the  reduced  frequency.  Since  the  aerodynamic 
derivatives  are  complex  quantities,  the  determinant  can  be 
expressed  into  two  parts,  real  and  imaginary.  After  expanding  the 
flutter  determinant  and  substituting  appropriate  values  of  Eqs. (42) 
the  real  and  imaginary  parts  of  the  flutter  determinant  become. 


" ° (44) 


f 


r 


* 


UK,Ko)p2  - UKZ(I  + B±1  CmaR)  + KQ(M£  - Mi  C4zR)]£ 
+ I (M£I  - S2  £)  - Ml 


0) 


2<lctzR-MllIcm<»R  + stCmzR  - stctaR> 


Pu»  (CizRCmnR  " cJtzICmaI  + ceaICmzI  " ctaRCmzR)1  ” 0 


[KaCPZI  - ‘’Vmd1  P + 'stc*,I  - ICtzI  + M'!c, 


nial 


(45) 


- S£C  T + 
mzl 


p£ 


to 


2 (C£aRCmzI  + C£aICmzR  “ C£zRCmaI  ” C£zICmaR) 1 0 


where  p = If  the  real  and  imaginary  parts  are  set  to  zero, 

P 

solutions  PI  and  P2  can  be  obtained  for  the  real  part  and  P3  for 
the  imaginary  part.  For  an  assumed  value  of  w , if  one  of  the 
solutions,  PI  or  P2,  of  the  real  part  is  equal  to  the  solution  P3 
of  the  imaginary  part,  then  this  to  corresponds  to  the  flutter  fre- 
quency. However,  the  aerodynamic  derivatives  are  functions  of  the 
Mach  number  and  the  reduced  frequency.  The  flutter  problem  can  only 
be  solved  after  the  aerodynamic  derivatives  are  evaluated.  It  is 
necessary  to  assume  a Mach  number  and  a reduced  frequency  and  test 
whether  flutter  occurs  at  these  values.  If  the  test  results  are 
negative,  then  it  is  necessary  to  iterate  on  reduced  frequency 
until  a flutter  case  occurs  at  these  values.  If  the  two  speeds 
are  not  equal,  then  it  is  necessary  to  iterate  on  Mach  number  until 
the  flutter  Mach  number  is  equal  to  the  assumed  Mach  number. 


RESULTS  AND  DISCUSSION 

A single  degree-of- freedom  analysis  for  torsion  was  performed  on 
cascade  of  blades  using  the  lifting  surface  aerodynamic  theory 
discussed  in  a previous  section.  After  obtaining  the  results, 
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a comparison  was  made  to  the  results  obtained  by  Whitehead  who 

12 

employed  a technique  developed  by  Smith  to  obtain  aerodynamic 
forces  and  moments.  Whitehead  distinguished  a region  of  flutter 
that  he  termed  as  "sub-critical  flutter",  which  occurs  in  a regime 
where  any  acoustic  waves  generated  cannot  propagate  upstream  and 
downstream.  Figure  1 shows  the  results  that  Whitehead  obtained 
for  a space- to-chord  ratio  of  1.0,  a stagger  angle  of  45°,  a position 
of  the  torsional  axis  at  58.8%  chord,  and  an  interblade  phase  angle 
of  60°.  The  imaginary  part  of  the  moment  coefficient  is  plotted 
against  the  reduced  frequency  parameter,  w . Also,  shown  in  Fig.l, 
is  a plot  of  the  imaginary  part  of  the  moment  coefficient  from  the 
present  calculations.  As  can  be  seen,  the  results  obtained  by  the 
Rao/Jones  lifting  surface  technique  agree  closely  with  the  reults 
obtained  by  Whitehead.  The  imaginary  part  of  the  moment  coefficient 
is  capable  of  adding  energy  to  the  system  when  it  is  negative. 
Assuming  that  there  is  no  mechanical  damping  involved,  flutter  will 
occur  as  a net  result  of  energy  being  added  to  the  system.  How- 
ever, in  reporting  his  results,  Whitehead  found  that  he  had  signifi- 
cant mechanical  damping  and  had  to  allow  for  an  average  value  as 
a flutter  limit.  This  limit  is  shown  in  Fig.l.  If  the  value  of  the 
imaginary  part  of  the  moment  coefficient  is  below  this  line,  the 
flutter  is  predicted.  The  point  where  it  is  just  possible  for  the 
flutter  to  exist  is  referred  to  as  a flutter  boundary.  Whitehead 
observed  that  as  the  Mach  number  is  increased,  the  points  at  which 
flutter  is  just  possible  move  to  progressively  lower  values  of  the 
frequency  parameter.  This  corresponds  to  higher  fluid  velocities 
or  to  lower  blade  stiffnesses.  The  effect  of  Mach  number  is  there- 
fore highly  significant.  i 

Figure  2 shows  the  frequency  parameter  below  which  torsional  flutter 
is  just  possible  and  has  been  plotted  against  blade  axis  position. 
Once  again,  a comparison  is  made  by  utilizing  the  Rao/Jones  technique 

and  Whitehead/ Smith  technique.  The  results  are  in  close  agreement 
and  concur  with  Whitehead's  observation  that  increasing  the  Mach 


number  is  favorable  and  that  the  worst  position  for  the  torsional 
axis  is  around  fifty  to  seventy  percent  chord. 

A numerical  investigation  for  two  degrees-of-f reedom  was  conducted 
using  the  unsteady  airlods  program  interfaced  with  the  flutter 
analysis  program.  The  geometric  and  structural  properties  of  the 
cascade  reference  airfoil  used  in  this  analysis  were  calculated 
from  the  data  obtained  from  Refs. 13,  14,  and  15.  The  properties 
used  are: 


Static  moment  about 
the  elastic  axis 

Moment  of  inertia  about 
the  elastic  axis 


S = 5.35xl0~4slug-ft 


I 

M 


6.316x10  ^slug-ft^ 


0)  = 


u>  = 


0.0214  slugs 
1000  rad/sec 
2000  rad/sec 
0.0833  ft 


Mass  per  unit  span 
Bending  natural  frequency 
Torsional  natural  frequency 
Blade  semi-chord  length 
Figure  3 represents  the  results  obtained  for  the  flutter  analysis 
for  various  values  of  interblade  spacing  (s  = 1.0*-,  1.1*-,  1.2*-,  and 
1.3£)  while  the  interblade  stagger  angle  was  varied  from  44°  to  60°. 
Over  the  range  of  values  considered,  it  can  be  seen  from  Fig. 3 
that  as  the  interblade  stagger  angle  is  i. .creased,  the  flutter 
Mach  number  also  increases.  This  indicates  that  an  increase  in 
the  interblade  angle  would  be  beneficial  in  preventing  flutter  at 
the  lower  Mach  numbers.  This  trend  is  in  general  agreement  with 
the  results  reported  in  Ref.  16. 

Figure  4 shows  the  effect  of  varying  the  interblade  spacing  for 
various  interblade  stagger  angles  (X  = 46°,  52°,  56°,  and  60°). 

The  trend  of  the  curves  shows  that  as  the  interblade  spacing 
increases  for  a given  interblade  stagger  angle,  the  flutter  Mach 
number  decreases.  This  indicates  that  increasing  the  interblade 
spacing  will  have  a destabilizing  effect  on  the  cascade  of  blades. 
This  same  trend  was  also  observed  in  Ref.  16. 

A summary  of  the  results  of  the  two  degree-of-freedom  flutter 
analysis  is  presented  in  the  Table. 


I Interblade  Spacing 


Figure  4 


Effect  of  Interblade  Spacing  on 
Flutter  Speed 


TABLE  Summary  of  Two  Degree-of-Freedom 
Flutter  Analysis 


Interblade 

Interblade 

Flutter 

Flutter 

Interblade 

Stagger 

Spacing 

Mach 

Reduced 

Phase 

Angle 

Number 

Frequency 

Lag 

X 

s 

M 

0) 

o 

(Degrees) 

(Degrees) 

1 . 21 
1.2  i 
1.21 
1. 3£ 
1. 3fc 
1.3A 
1 . 31 


0.4087 

0.5028 

0.5497 

0.5676 

0.6109 

0.6216 

0.6300 

0.6445 

0.3783 

0.5561 

0.5843 

0.5818 

0.6058 

0.4300 

0.4809 

0.5336 

0.5406 

0.5501 

0.4045 

0.4715 

0.4757 

0.4990 


0.4539 

0.3671 

0.3347 

0.3234 

0.2988 

0.2935 

0.2888 

0.2809 

0.4907 

0.3301 

0.3128 

0.3138 

0.3001 

0.4305 

0.3837 

0.3440 

0.3390 

0.3325 

0.4580 

0.3911 

0.3872 

0.3682 


CONCLUSIONS  AND  RECOMMENDATIONS 


A unique  velocity  potential  lifting  surface  technique  has  been 
developed  and  used  to  determine  the  unsteady  airloads  for  an 
infinite  cascade  of  staggered  blades  in  subsonic  flow.  The 
unsteady  airloads  were  then  utilized  in  a FORTRAN  program  to 
perform  a flutter  analysis  for  a single  degree-of- freedom  system 
in  torsion  and  for  a two  degree-of-freedom  system  in  bending  and 
torsion. 

Several  geometric  and  flow  parameters  of  a staggered  cascade  were 
varied  over  a specific  range  of  values  and  were  found  to  influence 
the  cascade  in  the  following  manner: 

1.  If  Mach  number  is  increased,  it  will  have  an  increased 
damping  effect  on  the  stability  of  the  cascade  in 
pitching  motion. 

2.  If  interblade  spacing  is  increased,  the  damping  decreases 
while  still  below  the  first  critical  frequency  is  equal 
to  zero  at  the  first  critical  frequency. 

3.  As  interblade  stagger  angle  increases,  the  lowest  cri- 
tical frequency  decreases. 

4.  As  interblade  phase  lag  is  increased  with  increasing 
reduced  frequency,  the  damping  will  decrease  rapidly  to 
zero  for  an  interblade  phase  lag  of  180  degrees. 

When  the  cascade  was  analyzed  for  a single  degree-of-freedom 
system  in  torsion,  results  were  found  to  compare  favorably  with 
the  results  obtained  by  Whitehead  for  flutter  boundaries.  A com- 
parison of  the  lifting  surface  technique  employed  by  Rao  and  Jones 
and  a different  technique  utilized  by  Whitehead  in  determining 
torsional  flutter,  were  shown  to  be  in  agreement  in  indicating  that 
increasing  Mach  number  is  a desirable  condition  from  the  stand- 
point of  flutter. 

The  developed  computer  program  is  general  and  is  capable  of 
predicting  flutter  analysis  of  a two  degree-of-freedom  system  of 
an  infinite  cascade  in  subsonic  flow.  It  was  concluded  that 
increasing  interblade  stagger  angle  would  be  beneficial  in  pre- 
venting flutter  at  low  Mach  numbers  while  an  increase  in  interblade 
spacing  would  produce  a destabilizing  influence  on  the  cascade  of 
blades. 
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While  a new  flutter  program  was  developed  to  determine  the  cascade 
flutter  speed  and  frequency,  only  a limited  range  of  values  were 
analyzed.  The  program  uses  very  efficient  techniques  in  computing 
unsteady  airloads,  the  flutter  frequency,  and  the  flutter  speeds. 
Further  investigation  should  be  conducted  over  a large  range  of 
cascade  parameter  values  for  different  geometric  configurations, 
before  one  can  attempt  to  make  general  conclusions. 
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APPENDIX  - PUBLISHED  RESEARCH  PAPERS  AND  REPORTS 


1.  Kronenberger,  Jr.  L. , "Flutter  Analysis  of  a Cascade  of 
Staggered  Blades  in  Subsonic  Flow,"  M.S.  Thesis,  Texas 
A&M  University,  December  1977. 

The  purpose  of  the  report  is  to  utilize  a numerical  lifting 
surface  theory  developed  by  Rao  and  Jones  to  predict  the 
unsteady  airloads  on  a cascade  of  staggered  blades  in  sub- 
sonic compressible  flow.  An  investigation  is  conducted  to 
determine  the  effect  on  the  unsteady  airloads  when  parame- 
ters such  as  reduced  frequency,  interblade  staggered  angle 
interblade  spacing,  and  interblade  spacing,  and  interblade 
phase  lag  are  varied  over  a specific  range  of  values.  Once 
the  unsteady  airloads  have  been  determined,  they  are  used 
to  perform  an  aeroelastic  analysis  of  the  staggered  cascade 
for  a single  degree  of  freedom  in  torsion  and  a two  degree- 
of-freedom  system  in  bending  and  torsion. 

Results  of  the  single  degree  of  freedom  analysis  yield  flutter 
boundaries.  These  are  compared  to  results  obtained  by 
Whitehead  who  used  a different  technique  for  calculating  tne 
unsteady  airloads  on  a finite  cascade.  A new  general  flutter 
program  is  developed  for  the  two  degree-freedom-system.  The 
airloads  are  used  as  forcing  functions  in  the  resulting  two 
Lagrangean  equations  of  motion  representing  the  bending  and 
torsional  degrees  of  freedom.  The  iterative  procedure  of  the 
flutter  program  yields  the  flutter  frequency  and  speed  of  the 
cascade  reference  airfoil  as  a function  of  the  cascade  para- 
meters. 

2.  White,  G.P.,  "Flutter  Analysis  of  a Cascade  of  Rotor  Blades," 
Presented  at  the  AIAA  13th  Annual  Meeting  and  Technical  Display, 
January  1977;  also  presented  at  the  24th  Annual  AIAA  South-West 
Student  Paper  Competion  (won  the  First  Place  Award  in  the  Under- 
graduate Division),  April  1976. 

A classical  two-dimensional,  bending- torsion  flutter  analysis 
of  a reference  airfoil  in  a cascade  of  infinite  blades  is  per- 
formed. The  unsteady  airloads  on  the  reference  airfoil  are 
predicted  using  a numerical  lifting  surface  theory.  Several 
cascade  and  flow  parameters  such  as  interblade  spacing,  stagger, 
phase  angle  between  blades,  Mach  number,  and  frequency  are 
investigated.  The  bending-torsion  flutter  speed  of  the  cascaded 
reference  airfoil  is  studied  as  a function  of  the  cascade  and 
flow  parameters  and  the  results  are  compared  with  that  of  an 
isolated  airfoil. 

3.  Jones,  W.P.  and  Moore,  J.A.,  "Aerodynamic  Theory  for  a Cascade 

of  Oscillating  Airfoils  in  Subsonic  Flow,"  AIAA  Journal,  Vol.  14, 
No.  5,  May  1976,  pp.  601-605. 

The  theory  is  developed  in  terms  of  the  velocity  potential 
rather  than  the  acceleration  potential,  and  a brief  outline  is 
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given  of  the  simple  numerical  technique  used.  The  effects 
of  varying  airfoil  spacing,  frequency,  Mach  number,  stagger 
angle,  and  phase  difference  between  adjacent  airfoils  are 
discussed.  Particular  attention  is  given  to  variations  in  the 
aerodynamic  damping  for  pure  vertical  translational  and  pitching 
motions.  It  is  shown  that  the  translational  damping  can  become 
zero  at  certain  discrete  frequencies  but  that  it  never  becomes 
negative.  The  pitching  moment  damping,  however,  can  become  nega- 
tive over  a wide  range  of  frequenicies  of  practical  interest. 

The  airfoils  are  assumed  to  be  zero  mean  incidence. 

4.  Rao,  B.M.  and  Jones,  W.P.,  "Unsteady  Airloads  on  a Cascade  of 
Staggered  Blades  in  Subsonic  Flow,"  Presented  at  the  46th  Meeting 
of  the  AGARD  Propulsion  and  Energetics  Panel,  Unsteady  Phenomena 
in  Turbcmachinery,  September  1975.  AGARD-CP-177 . 

The  Jones-Moore  numerical  lifting  surface  technique  is  applied 
to  the  theory  developed  by  Jones  to  predict  the  airloads  and 
moments  on  an  airfoil  of  a staggered  cascade  of  rotor  blades  in 
subsonic  flow.  Circumferential  distortion  due  to  inlet  flow 
conditions  is  expressed  as  an  interblade  phase  lag  and  both 
cases  of  oscillating  airfoils  and  oscillatory  inflow  were  consi- 
dered. Results  obtained  for  several  values  of  frequency,  stagger 
angle,  blade  spacing,  and  interblade  phase  lag. 

5.  Rao,  B.M.  and  Jones,  W.P.,  "Unsteady  Airloads  on  a Cascade  of 
Staggered  Blades  in  Incompressible  Flow,"  Proceeding  of  a 
Workshop  on  Unsteady  Flows  in  Jet  Engines,  Sponsored  by  the 
Project  Squid  (US  Navy),  AFOSR,  and  UARL,  November  1974. 

The  Jones-Moore  numerical  lifting  surface  technique  is  applied 
to  predict  the  airloads  in  oscillatory  flow  on  an  airfoil  of  a 
staggered  cascade  in  incompressible  flow.  Circumferential  dis- 
torsion  due  to  inlet  flow  conditions  is  expressed  as  a phase  lag 
between  blades  as  suggested  by  Schorr  and  Reddy.  Also,  the 
results  are  obtained  for  one  combination  of  stagger  angle  and 
blade  spacing  at  several  phase  lag  angles  and  frequencies. 
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